Rays, modes, wavefield structure and wavefield stability 
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Sound propagation is considered in range-independent environments and environments consisting 
of a range-independent background on which a weak range-dependent perturbation is superimposed. 
Recent work on propagation of both types of environment, involving both ray- and mode-based 
wavefield descriptions, have focused on the importance of q, a ray-based "stability parameter," and 
/3, a mode-based "waveguide invariant." It is shown that, when f5 is evaluated using asymptotic 
mode theory, (3 = a. Using both ray and mode concepts, known results relating to the manner by 
which a (or /3) controls both the unperturbed wavefield structure and the stability of the perturbed 
wavefield are briefly reviewed. 

PACS numbers: 43.30.Bp, 43.30.Cq, 43.30.Ft 



I. INTRODUCTION 

In recent workiiSiiii^ on sound propagation in the 
ray limit in weakly range-dependent ocean environ- 
ments it has been shown that both ray amplitude 
and phase (travel time) distributions are largely con- 
trolled by a property — described by the stability param- 
eter a — of the background sound speed profile. Other 
investigatorslQ-iiii2.i3.i4.i7.i8,i9,20^2i,22.23 j^aking use of 

a modal description of the wavefield in both range- 
independent environments and range-independent envi- 
ronments with weak perturbations superimposed have 
shown that many wavefield properties are controlled by 
a property — the waveguide invariant (3 — which is defined 
using mode-based quantities in the background environ- 
ment. In this letter it is shown that when (3 is evaluated 
using asymptotic mode-theoretical results, (3 — a. 

In the two sections that follow a and /3, respectively, 
are defined and briefly discussed. In the final section 
those properties of wavefields that are known to be con- 
trolled by a or are briefly reviewed. To keep our pre- 
sentation brief we show only equations required to derive 
the result /3 ~ a or those that provide insight into wave- 
field properties discussed in the final section. 



II. RAYS: a 

We consider underwater acoustic wavefields in three 
space dimensions that are excited by a point source. 
We assume that the environment consists of a depth- 
dependent background, with a weak perturbation su- 
perimposed, so the sound speed is written as c{z,r) — 
C{z) + 5c{z,r). Here —z is depth and the range r is the 
horizontal distance from the source. When (5c ^ it 
is assumed that azimuthal coupling of the wavefield is 
negligible, so that one need only consider propagation 



in the (r, z) plane. To make our discussion concrete we 
will assume that C{z) has a single minimum, but this 
assumption is not necessary. 

It is well known (cf. e.g. Refs. |^0) that substitution 
of the geometric ansatz 

?2(z,r,(T) = a(z,r)e''^'^(^'''), (1) 

where u{z,r,a) is the Fourier transform of the pres- 
sure u{z, r, t) and a is the acoustic frequency, into the 
Helmholtz equation and collecting terms in descending 
powers of a yields the eikonal and transport equations. 
The eikonal equation can be solved for the travel time T 
by integrating the ray equations, 



dp^ _ dH dz _ dH AT _ dz ^ 
dr dz ' dr dpz ' dr dr 



where 



H{pz,z,r) = -\Jc ^{z,r) 



Pi 



(2) 



(3) 



It follows from these equations and the relationship 
dz/dr — tant^, where ip is ray angle with respect to the 
horizontal, that cpz = sin (p, so the vertical slowness Pz 
can be thought of as a scaled angle variable. For a point 
source it is convenient to label rays by their pz value at 
{z,r) = (2:0,0), Pz,o- The Hamiltonian H = —pr, where 
Pz and pr are the vertical and horizontal components of 
the slowness vector p with ||p|| = c~^{z,r). In a range- 
independent environment c = C{z), pr is constant fol- 
lowing a ray. The transport equation can be reduced to 
a statement of the constancy of energy fiux in ray tubes; 
its solution, accounting for azimuthal spreading, can be 
written 



al^\H{dz/dpz,o)r\ 



(4) 



Here Og is the value of a" 



"■^ at the small distance (1 m by 
convention) tq from the source. The partial derivative in 
Q is evaluated keeping r fixed: z{pzfi,r) is ray depth. 
In a range-independent environment an alternative form 
of (gl) is 



^'Author to whom correspondence should be addressed. Electronic 
mail: mbrown@rsmas.miami.edu 



,2 ' 



\PrlPz 



r \pz (dR/dpr) 



(5) 



Typeset by REVTeX 



2 



where R{pr, z) is the range of a ray, and ray depth is held 
constant in the partial derivative. 

For the range- independent problem z(r) and Pz(r) (fol- 
lowing rays) are periodic functions. This periodic motion 
is most naturally described using action-angle variables 
{!,'&)■ The transformed ray equations (cf. e.g. Refs. 00 
for details of the transformation) are 



dl _ dH _^ 
dr dd ' 



(6) 

iLuil) - H{T) (8) 
where H{I) is the transformed Hamiltonian. The action 



di9 _ dH 

dr 
dT 

dr 



1 



/=— / dzp,{z)^- r dz^C-^{z)-pl, (9) 

where z± correspond to the ray upper (-I-) and lower (— ) 
turning depths where C~^{z-) = C"^(z+) = Pr- The 
angle variable d increases by 2-k each time a ray completes 
a cycle, and = 27r/i?^(J) where Ri is the range of a 
ray cycle (double loop). 



ReiPr) 



dl 

dpr 



2pr 



dz 



_ ^C-^{z)-pI 



(10) 



where I{pr) is defined in (|5J|. The action-angle form of 
the ray equations can trivially be integrated: 



'&{r) = 130 + ^(1^ mod27r 
I{r) = h 

T{r) = [Ml)-H{I)]r. 



(11) 
(12) 

(13) 



[More correctly, a term d{G — Ii}) /dr should be added 
to the r.h.s. of © where G is the generating function of 
the canonical transformation {pz,z) — > (/,??). The corre- 
sponding endpoint corrections to H13|l is both numerically 
insignificant and not relevant to the discussion that fol- 
lows.] 

The stability parameter a{I), whose relationship to 
wavefield properties will be reviewed in the final section, 
is defined as 



ad) 



I duj 

Lu{I)'dl' 



It follows from the relationship uj{I) 
a{I) can be expressed in the form 



aipr) — 27r 



IjPr) dRe 
Rj{I) dpr 



(14) 

2tt/Ri{I) that 
(15) 



III. MODES: /3 

In a stratified environment c = C{z) in three space 
dimensions the modal decomposition of the wavefield has 



the form 



u{z, r, a) = 



E 



{zo,a)(l)miz,a)H^ '{ ^Prm^ ) 
Jdz (/'^(z,cr) 



(16) 



Here Hq'^'' is the zeroth-order Hankel function of the first 
kind and the normal modes satisfy 



dz2 



(17) 



together with a pair of boundary conditions. Here 
{aprmf' is a separation constant. For sound speed pro- 
files with a single minimum it is well known (cf. e.g. 
Ref. i for the outline of a uniform asymptotic deriva- 
tion) that an asymptotic analysis of H17|l for modes with 
turning depths within the water column reveals that each 
<b^ is associated with a discrete value of the action /, 



0-I{Prrn) = m + i, TO = 0, 1, 2, ■ 



(18) 



where I{pr) is defined in 0. 

We now introduce, following the argument given 
in Ref. 9, the important concepts of modal 
group and phase slowness. Owing to the orthog- 
onality of the modes, the quantity Um{r,a) = 
J dzu{z,r,a)(j)jyi{z,a)/ J dz4>1^{z,a) isolates the contri- 
bution to the wavefield from the mode with frequency a 
and mode number to. The inverse Fourier transformation 
of UmiT, f), weighted by s(cr), the Fourier transformation 
of the source time history s(i), is denoted u„i{r, t). If s{a) 
has a narrow bandwidth, centered at do 7 then a Taylor 
series expansion of krra = oprm, with prm — Prm{o') via 
Eq. (fT^ . about (Tq yields the result 



(19) 



where the envelope function '^mij', t) travels at the group 
speed, (dkr/da)^^ , evaluated at the center frequency and 
mode number to. The group slowness is defined as 



_ dkr 



(20) 



Note that H19|l represents a slowly varying dispersive 
wavetrain whose envelope moves at the group slowness, 
but within which surfaces of constant phase move at the 
phase slowness kr/a = Pr- 

Consistent with the asymptotic analysis presented here 



SgiPr) 



ReiPr) 



(21) 



(cf. e.g. Ref. IHorli the latter reference also includes, 
with additional references, the exact expression for 5*^). 
Here Re{pr) is given in (fTn|l and Ti{pr) is the correspond- 
ing expression for the single-cycle travel time. 



Te{pr) = 27Tl{pr) +PrRl{Pr) 



= 2 



^c-Hz)-p^, 



(22) 
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Note that although, according to H21() . Sg depends only 
on the turning depths of a mode, this dictates via the 
quantization condition (|18|l that Sg is in general a func- 
tion of both frequency a and mode number m. 

The "waveguide invariant whose relationship to wave- 
field properties will be reviewed in the next section, is 
defined as 



dSg 

dpr 



(23) 



where pr is the modal phase slowness (often written as 
Sp). Like Sg, in general /3 is a function of both a and m. 
It follows from Eqs. 1)101 121|[^^ that f3 can be expressed 



as 



I{Pr) dRe 

Rj{I) dpr 



a{pr). 



(24) 



This is the main result of this letter. It should be em- 
phasized that our modal analysis is based on asymptotic 
results, so that we have only demonstrated the asymp- 
totic equivalence of (3 and a. 



IV. WAVEFIELD STRUCTURE AND 
STABILITY 

In this section we briefly review those properties of 
wavefield structure and wavefield stability (to a small 
range-dependent perturbation) that are known to be 
controlled by a or /3. Fig. ^ shows wavefield intensity 
\u{z,r,t)\ in a range-independent sound channel in the 
depth-time plane at a fixed range, r — 500 km, for 
waves excited by a transient compact source with cen- 
ter frequency /o — 75 Hz and bandwidth A/ « 30 
Hz where a = 27r/. This plot was produced by solving 
the Thomson-Chapman^^ parabolic equation. Use of a 
parabolic approximation to the Helmholtz equation in- 
troduces some minor distortion of the wavefield, but the 
connections described here between a, /3, wavefield struc- 
ture and wavefield stability hold whether a parabolic ap- 
proximation is introduced or not. 

Perhaps the simplest interpretation of a is that it is 
a measure of the rate at which small elements of the 
extended phase space {z,p, T) are deformed by the back- 
ground flow (with dz/dr treated as the z-coordinate of 
a fluid motion, etc.) by shearing motionjiSi Although 
this property does not correspond directly to any ob- 
servable wavefield feature, this property helps to im- 
derstand other observable wavefield features. A very 
simple observable wavefield property controlled by a is 
travel time dispersion. It follows from the third of Eqs. 
© that dT/d/ = /(dcj/d/)r = a{I)uj{I)r. In a range- 
independent environment / is a ray label that increases 
monotonically with increasing axial ray angle, so this 
equation describes the rate of change of ray travel time 
with increasing axial ray angle. This is illustrated in the 
ray simulation in Fig. ^ "where a < except in small an- 
gular bands corresponding to steep rays; note that this 




FIG. 1: Upper panel: Intensity of a transient acoustic wave- 
field in a range-independent deep-ocean environment as a 
function of depth and time at a fixed range of 500 km. The 
wavefield was excited by a compact source near the sound 
channel axis with /o = 75 Hz and A/ « 30 Hz. The dynamic 
range of intensity is 40 dB. Middle panel: Corresponding ray 
travel times vs. depth. The light (heavy) line corresponds to 
rays with an upward (downward) inclination at the source. 
Lower panel: Corresponding plots of a vs. T, computed by 
eliminating the /-dependence from a{I) and T{I;r). 



ray property is also clearly visible in the corresponding 
finite frequency wavefield. Note also that zeros of a cor- 
respond approximately to cusps in the (z,T) plane. At 
such cusps geometric amplitudes diverge. This can be 
seen from Eq. (O. In that expression R is the total 
range of a ray, which can be written as nRi (n complete 
ray loops) plus end-segment corrections which depend on 
source and receiver depths and ray inclination (positive 
or negative) at the source and receiver. At long range 
the dominant contribution to R is nRg. Then it follows 
from (O and H15|) that at caustics either pz = or a = 0. 
[We emphasize that is true only in the large r asymptotic 
limit. At short range {dR/dpr)z has complicated struc- 
ture, usually including a singularity when = 0, so that 
application of © generally requires great care near such 
points.] The same argument reveals that at long range 
geometric amplitudes are inversely proportional to jaj ; 
this property is also seen (away from caustics and inter- 
ference fringes) in the finite frequency wavefield shown in 

Fig.n 

In the presence of a weak range-dependent perturba- 
tion a has been shown to control both ray stability^ 
quantified by Lyapunov exponents, for instance, and sev- 
eral measures of travel time spreads^. The latter prop- 
erty is linked to the travel time dispersion property in 
the background environment noted above. Also, in the 
presence of a weak perturbation, a plays a critical role in 
ideas relating to ray dynamics 4*^*^*2^ References 11120 31141 
show that unperturbed rays for which \a\ is large are 



generally more sensitive to environmental perturbations 
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than those for which \a\ is small. 

In a range-independent environment (3 controls modal 
group delay dispersion: modal group delays statisfy 
Tg = SgV, where Tg — Tg{m,a) and Sg = Sg{m,a), 
so {dTgl dpr)r = —(3r. The property that is most com- 
monly usedi^ to motivate the introduction of the waveg- 
uide invariant /3 is that, following a surface of constant 
wavefield intensity, a relationship 5r/r = (35(7 /a must be 
satisfied. This leads naturally to the defining relation- 
ship (|23|l . (Note that in some references, including the 
original work by ChuprovjiS —dSg/dpr is defined to be 
f//3 rather than (3.) This defining relationship H23I) im- 
plies a power-law dependence of the difference in eigenval- 
ues of the Helmholtz equation on the acoustic frequency, 
km — krm = cimnO'~^ where m and n are mode numbers 
and amn is a constant. More generally Grachev^^ argued 
that km — krm = bmnil^^^ whcrc rj is a parameter, 
such as water depth, describing the environment and 7 
is a constant. Consistent with this dependence is the re- 
lationship Sr/r — (35u /a + {P/j)6r]/rj. It is important to 
note (cf. also Ref. 14) that because /3 = /3(cr,m), when 
modes corresponding to different values of /? are not tem- 
porally resolved)^ application of the above results may 
be difficult. Much work relating to P has focused on 
the homogeneous, constant depth waveguide for which (3 
is independent of a and m for modes corresponding to 
small values of axial ray angle; for that problem /? can 
be approximated as a constant. 

In an environment consisting of a range-indenepent 
background on which a range-dependent perturbation 
is superimposed, (3 controls the spread of modal group 
delays owing to mode coupling. This is seen by not- 
ing that the total group delay of modal energy that has 
been scattered among mode numbers m^, i — 1, • • • , N, 
is Tg = J2i Sg{mi, a)ri where r — is the range. It 



follows that 

Tg = Sg{fh,a)r+^{fh,(j)^^{h-I)r, 

= Sg{m,a)r + f3{m,a)-^y2{I, ~ I)n {25) 

where / = (m -I- \)/<J = I{pr) amd m is a suitably cho- 
sen mode number (e.g. the mode number at r = or the 
average mode number). Also in range-dependent envi- 
ronments, the relationship 5r/r = f35a/a -\- {P/j)d'i]/ri 
has been used in the interpretation of measurements 
in shallow water waveguides^SiSiiSSiS^ and time reversal 
applicationsii^iiiiiSiiS Note that large \f3\ is associated 
with a high degree of wavefield sensitivity, e.g. high sen- 
sitivity of modal group delays to an environmental per- 
turbation, or high sensitivity of the location of a time 
reversal focus to an environmental perturbation Sr]. 

In retrospect it is not surprising that there should be a 
simple connection between a and (3 inasmuch as, whether 
one adopts a ray or modal wavefield description, the ob- 
ject of study, the wavefield, is the same. The result 6 ~ a 
is an aspect of ray-mode duality (cf. e.g. Ref. Il5f) . We 
expect that there are many wavefield properties that we 
have overlooked that are controlled by this parameter. 
We are unaware of any parameter that characterizes an 
acoustic environment whose important rivals that of a 
(or /?). 
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